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ABSTRACT 

Observations suggest that accretion discs in many X-ray binaries are likely flared. An outer 
edge of the disc intercepts radiation from the central X-ray source. Part of that radiation is 
absorbed and reemitted in the optical/UV spectral ranges. However, a large fraction of that 
radiation is reflected and appears in the broad-band X-ray spectrum as a Compton reflection 
bump. This radiation is delayed and variability is somewhat smeared compared with the in- 
trinsic X-ray radiation. We compute response functions for flat and flared accretion discs and 
for isotropic and anisotropic X-ray sources. A simple approximation for the response function 
which is valid in the broad range of the disc shapes and inclinations, inner and outer radii, and 
the plasma bulk velocity is proposed. We also study the impact of the X-ray reprocessing on 
temporal characteristics of X-ray binaries such as the power spectral density, auto- and cross- 
correlation functions, and time/phase lags. We propose a reprocessing model which explains 
the secondary peaks in the phase lag Fourier spectra observed in Cyg X-1 and other Galac- 
tic black hole sources. The position of the peaks could be used to determine the size of the 
accretion disc. 
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1 INTRODUCTION 

The X-ray spectra of radio-quiet active galactic nuclei (AGNs) and 
X-ray binaries can often be decomposed into a power-law like con- 
tinuum with a cutoff at 100 keV (Zdziarski 1999), soft black body 
emission from the accretion disc and/or neutron star surface, and 
the Compton reflection continuum with the associated iron emis- 
sion line at ~ 6.4 keV (see e.g. George & Fabian 1991; Nandra 
& Pounds 1994; Reynolds 1999). The reflection features are pro- 
duced in a rather cold neutral material which is often identified 
with an accretion disc. Correlation between the amplitude of re- 
flection, R, and the photon spectral index, T, in X-ray binaries and 
AGNs (Zdziarski, Lubinski & Smith 1999; Gilfanov, Churazov & 
Revnivtsev 2000a) implies that a relatively large fraction of reflec- 
tion originates close to the X-ray emitting region. In the case of X- 
ray binaries, some fraction of reflected photons, however, can also 
come from a companion star (Basko, Titarchuk & Sunyaev 1974; 
Done & Zycki 1999; Vikhlinin 1999) or an outer part of the accre- 
tion disc. In the case of AGNs, the delayed reflection from a distant 
molecular torus (e.g. Ghisellini, Haardt & Matt 1994) can be clearly 
observed when the central X-ray source turns off (Guainazzi et al. 
1998). 

The reflected radiation is necessarily delayed relative to the 



E-mail: juri.poutanen@oulu.fi 



direct radiation from the X-ray source. Studies of these delays can 
help to determine the geometry of the accretion disc, the position 
of the X-ray source, and the distance to the companion (Vikhlinin 
1999). A number of papers were devoted to studies of the response 
of the Fe line profile to an X-ray flare in the vicinity of a black hole 
(e.g. Reynolds et al. 1999; Hartnoll & Blackman 2000). Observa- 
tions of AGNs with future X-ray missions such as Constellation- 
X may open a possibility of determining directly from the Fe line 
profile temporal evolution, the geometry of the X-ray emitting re- 
gion and the mass and spin of black holes (see Reynolds 1999; 
Fabian et al. 2000 for reviews). However, in order to do that, ob- 
servation of a bright individual flare is needed. Otherwise, the in- 
terpretation of the line profile produced by many flares would be 
complicated. Such studies are impossible in the case of X-ray bi- 
naries even with future instruments since the photon flux per light 
crossing time of one gravitational radius is thousands times smaller 
than in AGNs. As an alternative to studying the line profiles one 
can use the whole available statistics at different energy bands and 
analyse the response of the reflected continuum radiation to a vary- 
ing X-ray source. 

All the temporal studies mentioned above assumed that the 
accretion disc is flat. In contrast, a number of observations indi- 
cate that accretion discs are flared, i.e. geometrically thick at the 
outer edge (see Verbunt 1999 for a recent review). There observa- 
tions include obscurations of the central X-ray source (e.g. White & 
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Holt 1982) and the variability properties in the UV/optical/infrared 
spectral band in low-mass X-ray binaries (e.g. Mason & Cordova 
1982; Vrtilek et al. 1990), the delays between the optical/UV emis- 
sion and the X-rays observed in the X-ray bursters (Triimper et al. 
1985) and in the super-luminal black hole source GRO J1655-40 
(Hynes et al. 1998). The inferred height-to-radius ratio, H/vout, 
at the outer edge of the disc varies from 0.15 to 0.5. This is much 
larger than predicted by the standard accretion disc theory (Shakura 
& Sunyaev 1973). 

In this paper, we study the impact of the X-ray reflection from 
flat and flared accretion discs on temporal characteristics of the X- 
ray binaries and AGNs. In § ^, we present the general formalism 
for computation of the response functions and introduce approxi- 
mations that allow one to consider a linear reflection response. For- 
mulae describing the impact of reflection on the temporal character- 
istics of the composite (direct and reflected) signal are presented in 
§ ^. Response functions for the flat and flared discs and the result- 
ing phase lags are presented in § Q Comparison with observations 
and discussion is given in § |5| Conclusions are presented in § |^ 




Figure 1. Geometry of the problem. Note that vector o does not generally 
lie in the plane defined by other vectors. 



relative to the photons reaching the observer directly. Here p = 
d h and p = y^r^ + {z — h)^ is the distance between the source 
and the reflection point. 

The temporal evolution of the radiation reflected towards the 
observer can be described by the following relation: 



2 RESPONSE FROM THE ACCRETION DISC 
2.1 General formulation 

Let us consider an axially-symmetric accretion disc and assume for 
simplicity that the height of the disc surface above the equatorial 
plane is given by a power-law relation 



H{r/rout) 



(1) 



where H is the maximum disc height, rout is the disc radius, and 
Q > 1. The inner radius of the disc is Hn. Throughout this paper, 
distances will be measured in units of Rg = 2GM/(? and time- 
in units of Rg/c. For all the examples considered below, we take 
M = IOMq, but the results can be easily scaled to any mass. 

Let us introduce the Cartesian coordinate system with the ori- 
gin at the disc center. The x- and y-axes lie in the central disc plane, 
and the line of sight towards the observer is in the x — z plane in 
the direction o — (sin i, 0, cos i), where i is the inclination. The 
geometry of the problem is presented in Fig. |l|. An X-ray point 
source is situated at the disc axis at some height h from the disc 
plane, i.e. at the position h = (0, 0, h). Let us assume that the X- 
ray source angular distribution does not depend on time and photon 
energy. The direct (marked with a superscript D) monochromatic 
luminosity (per unit solid angle) is then 



Li(M,t)=Li(t)^. 

The total emitted luminosity is then 



(2) 



(3) 



Function gives the dependence of the incident (direct) 

radiation field on zenith angle arccos fi, and is normalised 
I J Q{n)dfi = 1, i.e. for an isotropic source, ^(/i) ~ 1. 

An element of the disc surface has coordinates d — 
(r cos (j), r sin (p, z), where the azimuth <j> is measured from the x- 
axis. Photons reflected from that element are delayed by 



L^^(cos i, t) = 



rdr I FE{ri,C,r,<t>,t)d(l>, 



(5) 



where Vm , ru are the innermost and outermost disc radii intersect- 
ing the isodelay paraboloid (see § 2.5 and Appendix A3) and 



FE{il,C,r,(t),t) 



»(m) C 2r? 
47rp^ cos ^ 27r 



j n{r-E,n;E'X)L''E,{t~M)dE' (6) 



is the monochromatic flux of radiation reflected in a unit solid angle 
from a surface element that gives a unit area when projected onto 
the central plane (note factor cos ^ = 1/ VT+z^ in the denomina- 
tor, where z' = dz /dr). The reflection Green function 71 describes 
the processes of Compton reflection, photoelectric absorption, and 
fluorescent line emission (see e.g. Poutanen, Nagendra & Svens- 
son 1996). It depends on the angles arccos arccos r) between the 
normal n to the disc surface at a reflection point and the directions 
of the incoming and outgoing photons, respectively, as well as on 
the photon energies. Here 



= — n ■ p/p = [r sin ^ — (jz — /i) cos ^]/p, 
77 = n ■ o = — sin j sin ^ cos + cos ^ cos j, 
11 = {z-h)/p. 



(7) 



At(r,( 



p — p ■ o 



r sm I cos ( 



h) cos i 



(4) 



In principle, there exists a weak dependence of TZ on the differ- 
ence in the azimuth of incoming and outgoing photons, which we 
neglect. This is equivalent to assuming isotropic scattering. The re- 
flection Green function can also be a function of radius, for exam- 
ple, if the ionisation state of the disc changes with radius. 

It is worth noticing that the evolution of the reflected radiation 
(eq. ^) depends not only on time, photon energy, and inclination 
angle, but also on the shape of the intrinsic spectrum L^{t). Thus, 
in general, temporal characteristics of the reflected radiation at a 
given photon energy E depend on behaviour of the intrinsic spec- 
trum at all energies above E. This means that there is no linear 
response between the intrinsic and reflected signals at a single en- 
ergy E. In a number of physically realistic situations considered 
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below the response is, however, very close to linear. That helps us 
to solve temporal problems in a simpler way. 

2.2 Constant intrinsic spectrum 

If the spectrum of the intrinsic X-ray radiation (t) does not vary 
in time while the normalisation varies, 

LiW = PWL°(i5), (8) 
equation can be simplified: 

Fein, C, r, t) ^L^it~ At)^^|^ps(r, r,, C), (9) 

Anp^ cos 4 Itt 

where we introduced the angle- and energy-dependent albedo func- 
tion 



PEir,ri,0= lnr;E,7T,E',0L°iE')dE'/L''{E). (10) 



The temporal and all other variables are separated. The temporal 
evolution of the reflected radiation at energy E and inclination i 
can be represented as a convolution of the direct luminosity at the 
same angle i with the energy- and inclination-dependent response 
(transfer) function: 



and is normalised to the reflection albedo: 



Lg"(cosi,t)= / TE{t-t')LE{cosi,t')dt' 



(11) 



Tsit) = rdr 5{t - At)d<t) 



cos ^ n(cos i) 2tv 



(12) 



The total observed luminosity is thus presented as a sum of the 
direct and reflected luminosities: 



PE{ri,02r]dri = asiO = 1 - Hx{OVT^. 



(16) 



For an isotropic source above the infinite plane the angle-averaged 
albedo then 



aE = j asiOdC, = 1 - /loVl - A, 



(17) 



where ho = 2(1 — yjl — A)/A is the zeroth moment of H\. The 
function H\ satisfies the integral equation 



+ (18) 



For A < 1, 

^AW^l + ^77lni±^, aE^). (19) 
2 r] 4 

One should note that the response function depends on the photon 
energy through the function pE- 



2.4 Isotropic approximation 

The response function in general depends on the properties of the 
Green function for reflection which is a function of the photon en- 
ergy and possibly radius (see eq. ^). However, if the physical con- 
ditions do not change with radius and if one assumes that the in- 
tensity of the reflected radiation is isotropic (we will call this the 
isotropic approximation), one can approximate pE{r, rj, C,) = aE- 
In this case, one can introduce a response function that is indepen- 
dent of albedo and photon energy: 



T(t).^= [ rdr [ 5it-At)d^^ 

^ ^ aE J J ^p^gog^ f^((,osj) 27r 



(20) 



Le (cos i,t) = L° (cos i,t) + / TE{t — t')L^ {cos i, t')dt' . (13) 



The advantage of this approximation is that we need to compute 
only one response function. In most applications considered in this 
paper we use this approximation. We discuss its accuracy in § 4.4. 1 



2.3 Thomson approximation 

X-ray observations with high temporal resolution and large photon 
statistics exist only at moderate energies (E ^ 20 keV), for which 
in the first approximation Compton recoil can be neglected, and the 
reflection Green function then takes the form 



n{r-EX\E' ,r^) ^ PE{r,7^,05{E ^ E'). 



(14) 



In the Thomson approximation (|l4|), the energy integral in equation 
(^ can be trivially taken. It is easy to see that the reflected radiation 
is then given by equation ( |ll] ) with the response function ([T^). 

A specific form of the function pE depends on the properties 
of the reflecting medium. If, for example, the medium is homoge- 
neous, i.e. the single-scattering albedo (the ratio of Thomson cross- 
section to the sum of Thomson and photoelectric cross-sections), 
A = (tt/ (fph + ctt), is constant with depth, the reflection Green 
function pE can be presented as a product of two Ambarzumian 
functions (see e.g. Chandrasekhar 1960; Sobolev 1975) 



Pe{vX) 



(15) 



2.5 Integration over azimuth and radius 

In order to integrate expressions ( [l^ and ( ^o| ) over the azimuth one 
can use the following identity: 



G{r,t) 



S{t- At)d(t> 



r sm I sm < 



(21) 



2/y/r'^ sin^i - [t - p + {z - h) cosi]2, i / Q, 
2TTS{t - p + z - h), i = 0. 

For given t and r, the azimuth is found from equation (^ (note 
that p and z are functions of r only) 



t — p + {z — h) cos i 



(22) 



The limits rm, rm in equation (|5|) are determined by the conditions 
that expression under square root in ( ^1] ) is positive (i.e. | cos 0| < 
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1, eq. g2|) and that rin < rm < '"m < ^out . The details on the 
method of integration over the radius in equation (||) are given in 
the Appendix. 



3 TEMPORAL CHARACTERISTICS 



As we showed in § 2.2 and § 2.3, the reflected signal depends lin 



early on the intrinsic (direct) signal when there is no spectral vari- 
ability of the direct radiation or at sufficiently low energies where 
Compton down-scattering can be ignored. Timing properties of the 
reflected radiation are then described by a response function Te it) . 
For most of our applications, we assume that albedo does not de- 
pend on radius and one can take TE{t) = aET{t), where now 
T{t) is energy independent. The total observed radiation flux is 
composed of the direct radiation from the X-ray source plus the 
radiation reflected from the disc (see eq. |l^ ): 



LE{t) = L'^it) + aE / T{t-t')L'^{t')dt' 



(23) 



The cross-correlation function (CCF) of the total signal, Ci2{t) = 
J Li(t')L2(t' + t)dt', can be easily found substituting equa- 
tion (^231): 



Ci2(t) = CE.(t) + aia2 / CY2{t')A'^{t-t')dt' 



+ 



C?2it') [a2Tit - t') + aiT{t' - t)] dt' , (24) 



where (7^2 (t) is the CCF of the direct signal and {t) = 
J T{t')T{t' + t)dt' is the auto-correlation of the response func- 
tion. The auto-correlation function (ACF) can be obtained from the 
above expression: 



AEit) 



+ / A'iit')A'^{t~t')dt' 



+ as / A'l:{t')[T{t' + t) + T{t' -t)]dt', 



(25) 



where is the ACF of the direct signal. 

The Fourier transform of (^) reads 



(26) 



The cross-spectrum of the light curves at two energies Ei,E2, 
Ci2{f) = [Li(/)]*L2(/), can be expressed asaproduct 



Ci2(/) = CUncUf) = \C?2if)\e 



(27) 



where C^2{f) = [L? if)]* if) is the cross-spectrum of the di- 
rect radiation and 



Cf2{f) = 1 + a^a2\f{f)\^ + aif *(/) + a2f{f). 



(28) 



The phase lags of the direct radiation ipu{f) and that added by 
reflection ipn{f) are defined as the phases of (712 (/) and (7i2(/)i 
respectively. The "reflection" phase lag. 



{a2-ai)QT{f) 



1 + (ai + a2)UT{f) + aia2|r(/)is 



(29) 



is positive (hard lags) when a2 > ffli. Here 5R and 9 denote real 
and imaginary part of a complex variable. For small ai and small 
02 / T{t)dt, the phase lag is approximately: 

Mf)~ ia2-ai)^f{f). (30) 
The power density spectra (PDS) are related in the following way: 

\LEif)\' = [1 + al\f{f)\' + 2aEm.f)] ■ OD 

In many recent publications, the PDS is often normalised to the rel- 
ative rms of the signal (Miyamoto et al. 1991; Nowak et al. 1999a). 
We denote such a PDS P{f) - In this normalisation, we get 



[l + ae/ T{t)dty 



(32) 



4 RESULTS 

4.1 Simple response functions 

For illustration, let us specify a simple exponential response func- 
tion: 



_Rexp[-(t-r)/r]/r, t > r, 

0, t<T. 



(33) 



The corresponding Fourier transform is T(f) = i?exp(ja;)/(l — 
ix), where x — 27r/r. As a starting point let us take Ei = 3 keV 
where for the neutral reflector of solar abundances the albedo ai ~ 
0.014. This is the smallest energy available for timing analysis, for 
example, at RXTE. For ue ^ ai , the phase lag at energy E relative 
to El is then 



ta.n<fiE{f) 



fleff {x COS X + sinx) 
1 + x^ + Ocff (cos X — X sin x) 



(34) 



where a^H = Rue is the effective albedo. The phase lag grows 
linearly with frequency as ifi^if) ~ 2aaax/{l + acff) for a; ^ 1, 
and for a^s Si 3 reaches the maximum tan <,Omax,i ~ O.Tflcff /(I — 
O.lSacff ) at a; ~ 1 (see Fig. ^). The peak shift slightly towards 
higher frequencies at larger tteff. When acff ^ 3 (large effective 
albedo is possible in strongly anisotropic sources), phase lag goes 
through TT at a; ~ 2. For ai <ti fleff ^ 1 and x <ti 1, the phase 
lags are proportional to the albedo ue (see eq. and their en- 
ergy dependence is exactly the same as that ob. At small energies 
close to El, lags go to zero. At larger energies (and larger albedos), 
the phase lag rises slower than aE (Fig. gb). At large frequencies, 
when cos a; ~ a; sin x (i.e. x ~ 1), the dependence on a_E in the de- 
nominator of eq. (Q) disappears and ip^ depends linearly on a^. 
Including the fluorescent iron line would produce a bump in the 
phase lag energy spectrum which amplitude depends on the energy 
resolution of the X-ray instrument. 

If the response function is a sum of two exponentials described 
by Ti, Ri and T2,R2, respectively (let ti ^ T2), the phase lag has 
two prominent maxima at /i ~ l/(27rri) and /2 ~ l/(27rT2). The 
value of the high frequency maximum at /i is the same as for the 
case of a single reflector, while the new, low frequency maximum 
is tan (^JJ;;,x,2 ~ 0.7acff,2/[l + a^ff,! - 0.15aoff,2]- 

Let us now consider the effect of reflection on the shape of 
the auto- and cross-correlation functions. The effect is easier to 
understand if one considers an even simpler 5-function response 
T{t) = RS{t - t). Then 

AE{t) = (l + a|i?')Ag(t)+ai5i?[Ag(t-r)+Ag(t + r)].(35) 
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Frequency, x=27t/t Energy [keV] 

Figure 2. Phase lags for the exponential response function (^sj). (a) 'fi^(f) frequency dependence for different a^ff relative to 3 keV (where ai = 0.014). 
Thin curves give the negative phase lags, (b) Solid curve gives the angle-averaged albedo for Compton reflection from neutral medium (omitting Fe Ka 
line) coiTesponding to the intrinsic photon power-law index F = 2 computed for solar abundances using pexrav model (Magdziarz & Zdziarski 1995) from 
XSPEC. Dotted curves give the energy dependence of the phase lag ^p^if) (relative to 3 keV) divided by frequency x = 27r/T for /? = 1 at frequencies 
a; = 0.1, 0.5, 0.7, 1, 1.5 (from top to bottom). Atx< 0.1, the curves coincide with the a; = 0.1 case. 



For monotonic, strongly peaked at zero lag ACF (such as that of 
Cyg X-1, see Maccarone, Coppi & Poutanen 2000) and ae 1, 
the largest relative increase in the total ACF achieved at lag ~ r is 
~ 1 -f a_B_R/A§(r) (for the normalised to unity at zero lag ACF). 
This means that the changes are more significant when the charac- 
teristic decay time-scale of the ACF is smaller than or comparable 
to the characteristic delay due to reflection. Similar behaviour is 
expected for the exponential response function (^3|). 

The CCF for ai = and 5-function response takes the form: 



Ci2(t) = CSW + a2i?CP2(i- 



(36) 



Since the CCFs of most GBHs and AGNs are strongly peaked 
at zero lag (see e.g. Nolan et al. 1981; Papadakis & Lawrence 1995; 
Smith & Liang 1999; Lee et al. 2000; Maccarone et al. 2000), 
the total observed CCF should become more asymmetric for larger 
albedos (i.e. at larger energies). The largest relative deviation from 
the intrinsic CCF should be observed at lags t ~ r. In the case of 
the intrinsic CCF (^{^(t) — exp[— (t/to)'^], the relative change at 
large lags, |t| 2> t, is the following: 



cm 



1 + a2Jiexp 



rty 



(37) 



4.2 Flat disc response 

4.2.1 Infinite slab 

The response function for a flat disc and an isotropic source can 
be derived analytically from equation ( |2c| ) in the isotropic approx- 
imation (i.e. p — as). In this case, z = 0, cos^ = l,rj = cosi, 
C, = h/p, and 



Tslab (i) 



2 cos i 
2tv 



r dr h 



G{r,t), 



(38) 



where = + and G{r,t) is given by equation (pi]). The 
variable p should satisfy the condition p_ < p < p+, where 



p±=q±s, q = 



t — hcoai amiy^t{t — 2/icosi) 



(39) 



Substituting p = q + s cos 6 into equation (BS) we get: 



Tslab(t) 



2 cos i / hdp 



27!" J cosi-y/ (p — p-){p+ — p) 



AO 



2h Txq 



2h 

TT J (g + SCOS&)2 TT (§2 — s2^3/2 




(40) 



h{t — hcosi) 

— 2C0Sl— -T f- ; —TTTTTy t> ZhCOSl. 

(t2+/i2_2t/iCOSi)3/2' 

The response function is normalised in the following way: 



TsiabWdf = 2cosi. (41) 



4.2.2 Disc with a central hole 

In hot flow models for the X-ray production in accreting black 
holes, the X-ray source is situated inside a cold accretion disc trun- 
cated at some radius (see e.g. Esin et al. 1998; Poutanen 1998). The 
response then can be approximately evaluated assuming a point 
source at the axis of the disc with rin 2> 1 (see also Gilfanov, 
Churazov & Revnivtsev 2000b). A more general expression for the 
response function, for a non-zero inner radius and a finite outer ra- 
dius Tout, can be written similarly to equation (|4^: 

PM 



2h 

TT 



dp 



de 



2h 

TT J {q + s cos ( 

9m 



(42) 
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Figure 3. (a) Response functions T(t) (in tlie isotropic approximation), their Fourier amplitudes, and the phase lags for the flat disc with different inner radius 
Tin and Tout = 10'' -Rg. Since asymptotically at large t in flat discs T(t) ~ t~^, the responses are multiplied by . The phase lags are computed assuming 
albedos ai = 0, 02 = 1. Solid and dotted curves correspond to cosi = 1 and 0.7, respectively. Thin curves give negative phase lags, (b) Dependence on the 
bulk velocity /9 = v / c. Solid, dotted, and dashed curves correspond to /? = — 0.2, 0, 0.3, respectively. Parameters are Tin = 0, rout = cos i = 0.7. 

The negative /3 means the bulk velocity directed towards the disc. Rg/c = 10"* s in all simulations. 



where 

tmin = max(2/i COS i, ^ r?^ -^h? -Vh cos i 
t max = \/ r-Qut + + /i COS i + rout sin i , 
and the Umits 



Tin Sm 1 1 



Pm = max 

PM = mill ( V^out + h'^, P+ ) 
9m,M = arccos[(pm,M - <7)/s]. 



(43) 



(44) 



For yj/ r?^ + h^ + hcosi + rin sin i < t < \/ r'^^yt + h'^ + h cos i — 
rout sini, the response function is given by equation (p^j). Fig. ^ 
shows the response function and other timing characteristics for 
different inner radii nn. Since in flat discs T{t) ~ at large t 
(see eq. |p^]), the responses are multiplied by t^. Due to the ab- 
sence of the reflector in the direct vicinity of the flare, reflection 
acts as a low pass filter removing high frequency signal from the 
reflected radiation. This causes the maximum of the phase lag to 
shift to lower frequencies. For large rin, the covering factor of the 
reflector is small and therefore the amplitude of the phase lag de- 



4.2.3 Anisotropic sources 

Let us consider anisotropic sources of radiation. When the X-ray 
emitting plasma has some bulk velocity (3 = v/c away from or 
towards the disc, the angular distribution of the radiation is 



1 



74(1 - /?^)3 



7 = l/y/T^, 



(45) 



if the emission is isotropic in the plasma rest frame (see e.g. Ry- 
bicki & Lightman 1979). This formula is valid for both the bolo- 
metric luminosity as well as for the monochromatic luminos ity for 



a photon spectral index F = 2. Similarly to derivations in § 4.2.1 



we then can analytically compute the infinite slab response function 
substituting fl{~h/p) to equation (ECl) : 



— (l-/3cosi)^ 

TV 



de 



p^l+f3h/p)3 







2^1^ .,3Y^ (fc + l)(fc + 2) , 

— (1-^cosi) > i-Ph) 

TV '—^ 1 



de 



(q + scos 



/i(l-/3cosi)=^^(fc + l)(fe + 2)(-/3/i)'= (46) 

fe=0 

where are the Legendre polynomials and we used relation 3.661 
from Gradshteyn & Ryzhik (1980) to obtain the last expression. 
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Figure 4. (a) Energy dependence of the reflection albedo for a slab with purely scattering skin of Thomson thickness r-r atop of a neutral slab, (b) Energy 
dependence of the phase lags. An isotropic X-ray source is elevated at = 10 above an infinite "ionised" slab. The phase lags are computed relative to Ex = 3 
keV at Fourier frequency / = l/(10/i) (i.e. 100 Hz for h = lOiJg = 3 10^ cm). 




1 10 100 1000 10 ^ 10^ 10^ 10^ 10*^ 

Time Delay [ms] Frequency [Hz] 



Figure 5. (a) Response functions Te (t) for a flat slab with ionised skin which optical depth decreases with radius as tt (r) = ro{l + r/10/t)~3/2, where 
TO = 1, and h = 10 are assumed (here cos i = 0.7 and Rg = 3 10^ cm). Thick curves from the bottom to the top correspond to 3, 6, 12, 24 keV. Thin solid 
curve represents the response function with unity albedo over the whole suii'ace. (b) Corresponding phase lags relative to 3 keV. Thin solid curve represents 
the case of constant albedos a2 = 1, ai = 0. 



The responses together with the related phase lags are presented 
in Fig. Bp for a few /3. With increasing fJ, the response amplitude 
decreases since radiation is becoming more beamed away from the 
disc towards the observer. At small lags the effect is even more 
dramatic. This causes also a significant change in the amplitude of 
the phase lags at high frequencies. 



4.2.4 Ionised discs 

The intense incident X-ray flux can affect the ionisation state of 
the reflecting material changing its albedo. The surface layer of the 
accretion disc is heated up to the Compton temperature of radia- 
tion and is almost completely ionised (e.g. Nayakshin, Kazanas & 
Kallman 2000). We can consider such "skin" as a purely scatter- 



ing medium. A transition layer from the hot layer to almost neutral 
matter is very thin. Therefore, a two-phase model, with a hot layer 
atop of neutral material, can be adopted. We neglect smearing of 
the reflection features by thermal motions in the skin and consider 
pure Thomson scattering. 

One should point out here that in the isotropic approximation, 
the response functions T{t) given by equation (EQ) are not affected 
by ionisation (while the reflected luminosity is affected). However, 
the phase lags at given energies will change since they depend on 
the reflection albedo. If reflection from the underlying neutral ma- 
terial follows Lambert's law (i.e isotropic intensity), the total re- 
flection albedo can be expressed as follows (Sobolev 1975): 



2(1 



1 + |C+(1-|C) exp(-rT/C) 
4 + 3rT(l-a) 



(47) 
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where rx is Thomson optical thickness of the ionised layer, and a 
is the albedo for a neutral slab. The angle averaged albedo is then 



o(tt) 



(1-a 
(1-a 



7-3 exp(-rT) + (4 + 2.tt)E2{tt) 

2[4 + 3rT(l-a)] 
7 + exp(-10rT) 



2[4 + 3tt(1- 



(48) 



where E2 is the exponential integral of the 2nd order. A relative 
error of the approximation in the nominator is smaller than 1% for 
any tt. The energy dependence of the albedo for different tt is 
shown in Fig. ^a. We see that at lower energies the albedo is high, 
amin(TT) ~ 3tt/(4 + 3tt), Compared with the neutral material 
albedo, and therefore one cannot assume ai <C 1 when comput- 
ing the phase lags. The important quantity here is, however, not the 
albedo itself, but the difference between albedos at different ener- 
gies (see eqs. E9l BQ): 



a2(rT) — ai(rT) 



(aa - ai){l - [1 - exp(-10rT)]/8} 
[1 + 37t(1 - ai)/4][l + 37t(1 - a2)/4] 
a2 — a\ 



(1 + 37t/4)2 



(49) 



The first approximate formulae is accurate within 1%. The accuracy 
of the second one is better than 13% for tt < 0.5, and is better 
than 20% upto albedo 02 ~ 0.7 (i.e. S ~ 20 keV) for tt ~ 1. 
The resulting phase lags (relative to 3 keV) at / = l/(10/i) as a 
function of photon energy are shown in Fig. kp. They have very 
similar behaviour at other frequencies (see Fig. ^) and follow the 
02(tt) — ai(rT) dependence (see eq. p9| ). 

In real physical situations, the optical depth of the ionised skin 
could be a function of the radius. For a source of ionising radiation 
at the disc axis at height h, we can approximate the radial distribu- 
tion asrT(r) = ro(l+r/10/i)~^/^ (see e.g. Nay akshin 2000). The 
resulting response functions Tsit) computed from equation ( p^ ) in 
isotropic approximation are shown in Fig. Close to the center 
(small delays), the responses for different energies are close to each 
other, due to the fact that the effective albedo is large for all energies 
(see eq. |p^). At large delays, a strong energy dependence of the 
effective albedo results in large differences in the responses. The 
response is strongly suppressed at small energies, while at higher 
energies the shape of the response function is similar to the constant 
albedo case. The resulting phase lags (see Fig. ^p) are then smaller 
(comparing to the constant albedo case) at higher frequencies since 
the difference of the albedos becomes smaller (see eqs. ^ 
^^). Phase lags at low frequencies are hardly affected by ionisa- 
tion. 



4.3 Flared disc response 

An important parameter influencing the response of a flared disc is 
the a-parameter (see eq. |^ which describes the disc shape. In the 
general case, it is difficult to get useful analytical formulae for the 
disc response, therefore we compute the responses numerically (see 
eq. [ po| | and Appendix). We first consider an isotropic source above 
the disc and vary the a-parameter. The results are shown in Fig. |^. 
For ct = 1 the response does not differ much from that of the flat 
disc. For any a > 1, there appears an increased response from the 
outer part of the disc and the response function flattens out at large 
lags. Comparing with the flat discs (and flared discs with a — 1), 
the additional phase lags corresponding to the light crossing time 
to the outer part of the disc r ~ rout /c appears at low frequencies 
/ ~ l/(27rr). The larger the curvature (i.e. larger a), the more 



prominent the bump in the Fourier-frequency-dependent phase lag 
spectrum becomes. 

Changes of the accretion disc size affect mostly the lags at 
low frequencies corresponding to the response from the outer edge 
(Fig. ^). For larger rout, the secondary phase lag peak shifts to- 
wards lower frequencies with a constant amplitude. 

Variations in the plasma bulk velocity and corresponding 
changes in the angular distribution of the intrinsic radiation af- 
fect markedly the flared disc response functions and the phase lags. 
They are presented in Fig. ^ for a few /3. The most prominent ef- 
fect is the dramatic change in the amplitude of the response at small 
lags, i.e. from vicinity of the X-ray source. The response from the 
outer disc edge is affected less since for 77 ~ and mildly relativis- 
tic velocities, beaming is not so important. The resulting phase lags 
are affected in a similar way: the amplitude of the phase lags at high 
frequencies changes dramatically, while the lags at low frequencies 
related to the outer disc edge remain basically the same. 



4.4 Approximations 

4.4.1 Accuracy of isotropic approximation 

Let us now estimate the accuracy of the isotropic approximation. 
Using the exact expressions for pe-function (eq. |^ we computed 
the response functions for the two cases A = 0.1 and 1 correspond- 
ing to the photoelectric absorption and the scattering dominated 
regimes. The results are shown in Fig. ^. The difference is coming 
from different angular dependences of the reflection Green func- 
tion ([1^ at different A. 

For small A, H\ ~ 1 and p_b(cos i, (^)/a_B ~ l/(cosi + ("). 
At rather small inclinations, cos i = 0.7 — 1, the reflected intensity 
is reduced comparing to the isotropic case (pe/ole = 1) at small 
lags (photons coming underneath the flare) since ~ 1. At large 
delays, ~ and the reflected intensity is closer to the isotropic 
case. For small inclinations, cos z ~ 0, the reflection is significantly 
larger than that given by the isotropic approximation at large delays 
(i.e. small Q (see lower panel in Fig.|7). 

For A ~ 1, Hx{rj) ~ 1 + 2?; and p_b(cos i, ~ 
(1 + 2 cos i)(l + 2(^)/4(cos i + At small inclinations, the isotro- 
pic approximation underestimates (overestimates) slightly the re- 
sponse at small (large) lags. At large inclinations and large delays, 
the response exceeds that in the isotropic approximation. 

Thus, the isotropic approximation overestimates (underesti- 
mates) response at small (large) lags for most inclinations and for 
A ^ 1. For A ~ 1, the isotropic approximation is very accurate at 
small lags. At large lags, it overestimates (underestimates) slightly 
response at small (large) inclinations. 

The phase lags expected for the flat disc and computed using 
the accurate response function as well as in the isotropic approxi- 
mation are presented on the right panels in Fig. For A <C 1, the 
"isotropic" phase lags are larger for small inclinations and smaller 
for large inclinations than the accurate ones. For A ~ 1, the phase 
lags computed in the isotropic approximation are almost indistin- 
guishable from the exact ones (compare dashed and solid curves 
peaking at higher frequencies). 



4.4.2 Approximate response function 

Calculations of the response function from a flared disc are not 
entirely trivial. However, for many applications they are not even 
needed. One can approximate the exact results by much simpler 
response functions that have almost identical temporal properties. 
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Figure 6. Response functions, their Fourier amplitudes, and tlie phase lags for flared discs in the isotropic approximation. The phase lags are computed 
assuming albedos ai = 0, 02 = 1. Thin curves give negative phase lags, (a) Dependence on a. Solid, dotted, and dashed curves correspond to a = 1, 1.5, 3, 
respectively, (b) Dependence on the outer disc radius rout- Solid, dotted, and dashed curves correspond to rout = 10*, lO" '^, lO^Rg, respectively, (c) 
Dependence on the bulk velocity /B = v/c. Solid, dotted, and dashed curves correspond to (3 = —0.2, 0, 0.3, respectively. Everywhere cos i = 0.7, and in 
(a), (c) rout = lO^Rg. 



The reflection response from the vicinity of the X-ray source can 
be approximated by the flat disc response (|4c|), while the response 
from the outer edge can be approximated by a power-law: 



(50) 



where Tout (t) = c ^"^+^■'^(""^''^^ Wn <t < rout, and the con- 
stant c is found from the normalization condition J Tout {t)dt — 1. 
We find that the resulting phase lags are similar to the exact ones 
(for Tin 7^ 0) if one takes 



^min — 



V ^hi + + ^ cos i, 

(1 — x)2hcosi + x{^Jrf^ + /i^ + hcoai 



x>l, 
X <1, 



that from the outer edge, depend on the inclination, bulk velocity 
/3, and solid angle occupied by the outer disc edge expressed in 

H/Vout- 



where x = rin/^tani. The normalisation factors, i.e. the apparent 
amplitude of reflection from the disc area close to the source and 



Ri 



R2 



(51) 



= (1 



-fJcosi) 



H 

rout 
H 

Tout 



H/rout 



(1 — /3 cos i) 



r2(cos i) 
3 1 - (/3/2)ff/rout 
(1 - l3H/routy ■ 



(52) 



(53) 



A more accurate approximation can be introduced if one replaces 
the first term in equation ( ^ ) by the "anisotropic" response (|46|): 
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Time Delay [ms] Frequency [Hz] 

Figure 7. Response functions and phase lags for the flat disc of radius rout = 10^ with an isotropic source at elevation h = 10 computed for three 
inclinations cos i = 1, 0.7, 0.2 (from top to bottom). Dotted and dashed curves represent the response function TE(t) given by equation ( [l^ computed using 
the exact angle-dependent expression for pE froin equation ( [l5| ) for A = 0.1 and 1, respectively, and divided by the angle-averaged albedo for these A (i.e. by 
= 0.026 and 1, see eq. ||l7||). The coiresponding phase lags (computed assuining ai = 0) are also divided by a^- Thin curves give negative lags. Solid 
curves correspond to the response function (ECl) computed assuming isotropic intensity of the reflected radiation, i.e. = a^- The phase lags are computed 
taking the same angle-averaged albedos, aE = 0.026 and 1. For larger albedo, the phase lags peak shifts to higher frequency. 




10^ 10^ 10^ 10^ 10 ^ 10 ^ 10° 10^ 10^ 10^ 

Time Delay [ms] Frequency [Hz] 

Figure 8. (a) Flared disc response functions for a = 3, H/rout = 0.2, h = 10, cosi = 0.7, Tin = 10, and rout = lO**. Dotted curves give the exact 
responses for /3 = —0.2 (upper curve) and /3 = 0.2 (lower curve). Solid curves represent the approximate response function (^^. (b) Coiresponding phase 
lags (coinputed for ai = 0, 02 = 1). 
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T,pp,(t) = Tp{t) + R2Tout{t). (54) 

This approximation is compared with the accurate response func- 
tions for flared discs at Fig. ||. One sees that it reproduces rather 
well the response functions. The approximation slightly overesti- 
mate (by 30% - 50%) the position of the peak in the phase lag 
spectra resulting in a similar error in estimation of the disc radius. 



5 COMPARISON WITH OBSERVATIONS AND 
DISCUSSION 

The phase lag Fourier spectra of a number of GBHs show breaks 
and secondary peaks at low frequencies (Nowak et al. 1999a,b; 
Grove et al. 1998; see Fig. ^). The appearance of a secondary peak 
could be related to the variability properties of the intrinsic sig- 
nal. The peaks in the phase lag spectra approximately correspond 
to the peaks in the power density spectra in the fP(f) represen- 
tation (see e.g. Nowak et al. 1999a; Kotov, Churazov & Gilfanov 
2001) which can be produced by two independent processes. The 
secondary phase lag peaks can also be a result of reflection from 
the outer part of the disc (Poutanen 2001). The position of these 
peaks would then depend on the disc size (see Fig. ^). For a linear 
response, the energy dependence of the lags at a given frequency 
should have features characteristic of reflection: an Fe line, an edge, 
and hardening above 10 keV. The lags observed in Cyg X-1 at fre- 
quencies above 0.6 Hz do not show these features (Kotov et al. 
2001). In fact, the reflection features appear to have negative con- 
tribution ("anti-lags"), i.e. the lags are reduced at energies where 
reflection is expected to contribute. This could be interpreted as 
a non-linear, negative response. At lower frequencies (below 0.6 
Hz), the energy dependence of lags is not known. It is possible that 
there the reflection response is still positive and the lag energy de- 
pendence shows reflection features. 

Let us develop a toy model that is capable of describing the 
observed phase lags (and the PDS). We assume that the lags at 
high frequencies are intrinsic, i.e. related to the direct radiation. 
We describe the direct emission in terms of a simple modified shot 
noise model. The shot profile is taken to be a rising exponential 
exp{t/tE), t < 0. For a power law radiation spectrum with time- 
varying spectral index, the shot time-scale depends logarithmically 
on energy ts = to[l — b\n{E/Eo)] (see Poutanen & Fabian 1999; 
Kotov et al. 2001). Spectral evolution can for example be produced 
by evolution of magnetic flares (Poutanen & Fabian 1999). We con- 
sider b as a free parameter since its value is determined by the de- 
tails of spectral evolution which are unknown. The decreasing of 
ts with energy ensures that the ACF width decreases with E as 
observed (Maccarone et al. 2000). The resulting CCFs peak at zero 
lag since the shots at different energies are correlated and peak at 
the same time. For a given to, the cross-spectrum of the signal at 
two energies is C{f, to) oc tit2/[(l — i27r/ti)(l + i2Trft2)]. If to 
is distributed according to P{to), the total cross-spectrum becomes 
C{f) = J C'(/, to) f'(io)dio- We assume a power-law distribution 
P{to) oc tg^ between tmin and tmax- The lags observed from Cyg 
X-1 on 1996, October 22 (Nowak et al. 1999a) at high frequen- 
cies can be described by this model (dashed curves in Fig. ^ with 
tmax ~ 0.05 s,p= 1.5, b — 0.22, and any sufficiently small tmin. 
At / < l/(27rtniax), phase lags are proportional to /. With these 
parameters, the PDS of Cyg X-1 above 1 Hz is also well reproduces 
(p = 1.5 gives the slope of the PDS also equal to 1.5, see Pouta- 
nen & Fabian 1999). The PDS shape at lower frequencies can be 
reproduced for example by the flare avalanches. 



Let us now consider the effect of reflection by applying an 
approximate response function for a flared disc given by equa- 
tion p^l). A rather good fit to the data can be achieved for the disc 
inner and outer radii 30i?g and lOOOiig, respectively, and the flar- 
ing parameter H/rout ~ 0.3 (see Fig. ^ we assumed (3 = 0, 
cosi = 0.7, and a = 3). One can see that reflection, however, 
cannot reproduce the observed sharp features. For the given pa- 
rameters, the amplitude of reflection is close to that observed in 
Cyg X-1 (e.g. Gierlihski et al. 1997; Poutanen 1998). In order not 
to overproduce lags at high frequencies, the reflection should be 
suppressed at small distances from the central source. This can be 
achieved by truncation of the disc at ~ 30-Rg(as here), by ionisa- 
tion of the disc, or beaming of the radiation away from the disc. 
A rather small size of the disc is consistent with the fact that Cyg 
X-1 is probably accreting via wind from the high-mass compan- 
ion. The angular momentum of accreting material is then small and 
the disc is small (Illarionov & Sunyaev 1975; Beloborodov & Illari- 
onov 2001). In low-mass X-ray binaries where material accretes via 
Roche lobe overflow, the angular momentum is high and the disc 
is large. In such systems the phase lag secondary maxima should 
appear at lower frequencies as seems to be the case (Nowak et al. 
1999b; Grove et al. 1998). Thus, if the model is correct, the disc 
size can be determined from the position of the phase lag peak. The 
lags energy dependence at low frequencies can resolve the issue 
whether the phase lag breaks are due to reflection. 

Revnivtsev, Gilfanov & Churazov (1999) showed that the 
frequency-resolved spectra (proportional to the amplitude of the 
Fourier transforms L^{f) at a given frequency / as a function of 
E) of Cyg X-1 depends on frequency. At / > 20 Hz, the spec- 
tra are rather hard and contain no reflection features. At low fre- 
quencies (below 0.1 Hz), the spectra are much softer and show 
strong reflection features. If the response is linear, reflection modi- 
fies the frequency resolved spectra in the following way | (/) 1 ~ 
[1 + aE^f{f)] (see also eqs. |T]). The straightforward 
explanation of the observed effects as solely caused by the finite 
light crossing time (studied in the present paper) fails, since the 
delays due to reflection cannot change the slope of the underlying 
spectra. The observed behaviour implies that the intrinsic (Comp- 
tonized) frequency-resolved spectra, -£§(/), depend on frequency. 
In that case, the signal at different frequencies possibly is produced 
in physically different locations. For example, the high frequency 
signal can be produced closer to the black hole, while the low fre- 
quency signal-further away. This arrangement would be natural if 
the dominant frequencies are the Keplerian ones. However, phys- 
ically different locations also imply different reflection response 
functions thus complicating the situation. 

The X-ray reprocessing in the cool disc produces both the 
Compton reflection continuum with an Fe line and soft radiation. 
The reprocessed radiation from the outer disc emerges in the op- 
ticalAJV bands. The response function of this soft radiation to the 
variability of the intrinsic X-ray radiation would be similar to that 
of the Compton reflection continuum unless ionisation state of the 
reflecting material is a function of radius. Thus, the accretion disc 
response function obtained from the cross-correlation analysis of 
the optical and the X-ray light curves (see e.g. Hynes et al. 1998; 
O'Brien & Horne 2001) should also satisfy the constraints coming 
from the X-ray timing e.g. from the Fourier-frequency-dependent 
phase lags. If the response is known from the optical/X-ray data, 
then one can try to subtract the lags related to the reflection from 
the observed time lags in order to obtain the lags corresponding to 
the intrinsic (direct) signal only. This will not be easy, however, if 
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Figure 9. Phase lags between signals in the 8.2-14 keV and the 3.9-6 keV bands vs the 2-3.9 keV band observed vi'ilhRXTE 'm Cyg X-1 on 1996, October 22 
(Nowak et al. 1999a; Pottschmidt et al. 2000). Dotted curves correspond to the intrinsic lags described by a modified shot noise model (Poutanen & Fabian 
1999) with exponentially rising shots of time-scale to distributed between 1 ms and 0.05 s according to a power law P{to) oc with p = 1.5. At a given 
energy E the shot time-scale is Ie = to[l — 61n(E/Eo)]. where Eq = Z keV and 6 = 0.22. Since decreases with energy, the lags are hard. Dashed 
curves are the lags produced by reflection from a flared disc with the following parameters = 10, 7*out — 

10^, /i = 10 (in units flg = 3 10^ cm), 
H/rout = 0.3, a = 3, cos i = 0.7, and /3 = 0. Solid curves give the sum of the lags from the intrinsic signal and reflection. 



the reflection response is non-linear or the disc is ionised (then the 
UV and Compton reflection responses are different). 

Real accretion discs are not necessarily axisymmetric, but can 
be, for example, radiatively warped (e.g. Pringle 1996; Wijers & 
Pringle 1999) or warped due to the Bardeen-Peterson (1975) ef- 
fect. For non- axisymmetric discs the response function depends not 
only on the inclination, but also on the azimuthal angle. Blackman 
(1999) and HartnoU & Blackman (2000) studied the Fe Une profiles 
produced in flared axisymmetric and warped discs. It would be in- 
teresting to study temporal response of the iron line profile for such 
discs. If reflection is responsible for shaping the phase lags, they 
would change periodically depending on the azimuthal angle. 

We focused on the effects of Compton reflection on temporal 
characteristics and not much attention was paid to the iron emission 
line. With future instruments of high spectral resolution and large 
effective area, it will be possible to study the phase lags across the 
line profile. By combining the line profile with the energy depen- 
dence of the phase lags it will be possible to put better constraints 
on the accretion disc geometry in X-ray binaries. The narrow cores, 
produced in the material far away from the X-ray source, would 
be associated with large time delays. The broad wings should re- 
spond faster to the variations of the continuum and the lags should 
be smaller. However, there could be serious complications. We as- 
sumed everywhere that there is a linear response between variations 
of the continuum and the reflected radiation. In other words, prop- 
erties of the reflecting material are assumed to be constant indepen- 
dently of the X-ray flux. This assumption could be wrong since the 
ionisation state of the reflecting material depends on the illuminat- 
ing flux (see e.g. Rejmolds 2000; Nayakshin & Kazanas 2001). 

Another complication is that the angular distribution of the 



intrinsic radiation could change. This would be likely for exam- 
ple in the model with mildly-relativistic outflows (Beloborodov 
1999; Malzac, Beloborodov & Poutanen 2001). If the velocity of 
the plasma ejection is correlated with the energy dissipation rate in 
magnetic flares, one could expect that a large emitted luminosity 
corresponds to a larger ejection velocity, and more beaming away 
from the disc. This would decrease the amplitude of reflection at 
large fluxes. The resulting reflection response could thus be non- 
linear possibly reproducing negative lags introduced by reflection 
(Kotov et al. 2001). Further studies in that direction are in progress. 



6 CONCLUSIONS 

In the present paper, we have developed the formalism for com- 
puting the response functions for isotropic and anisotropic X-ray 
sources above flat and flared accretion discs. We have also stud- 
ied the impact of the reflection on the temporal characteristics such 
as the power density spectra, auto- and cross-correlation functions, 
and the Fourier-frequency-dependent time/phase lags. A number 
of useful approximations was introduced to simplify significantly 
computations of the response function. Simple approximate re- 
sponse function was shown to have temporal characteristics which 
are very similar to that computed using exact responses. 

The temporal characteristics predicted by a model with a lin- 
ear response of the reflected component onto the variabiUty of the 
intrinsic radiation were computed. The reflection model can re- 
produce the observed secondary maxima at low frequencies in the 
phase lag Fourier spectra of Cyg X-1 and other GBHs. The position 
of these maxima could be used to get constraints on the accretion 
disc size. In spite of the fact that the observed energy dependence 
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of the phase lags at high frequencies indicate that reflection is not 
their only source (Kotov et al. 2001), a proper account of the im- 
pact of the reflection on temporal characteristics is necessary, when 
interpreting the data. 
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APPENDIX A I 
FUNCTION 



CALCULATION OF THE RESPONSE 



Al Visibility conditions 

Contrary to the case of the flat disc, the every point of a flared disc 
is not necessarily visible. For a concave disc surface, z" > 0, the 
condition for visibility reads 



tani < 



\/ ''out ~ sin'^ (j) ^ cos 4> 



H 



For a given inclination i and radius r, this translates to 



cos m < cos (ho 



(r^ 

\' OU 



2{H — z)r tan i 



(Al) 



(A2) 



According to the value of tani, the visibility condition is: (i) 
if tani < l/^max. then the whole disc surface is visible; (ii) 
if l/^max < tani < Vout/H, then the whole area r < Vc 
is visible. For rc < r < rout, the condition cos(j} < cos 00 
should be fulfilled. Here rc is the solution of the equation tan i — 
(rout — r)/{H — z); (iii) if tani > Vout/H, then the area r < rc 
is invisible. For r > rc, the inequality cos(/) < cos 00 gives 
the visible area. Here rc is the solution of the equation tani — 
(rout +r)/(i/ - z). 
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Figure 10. Upper panels: The curves of equal delays (isochrons) at the disc surface for different inclinations (only one half of the disc is shown). The 
reflector's surface is given by equation (|lj) with H/rout = 0.3, a = 3, rout = 1000, = 0. The flare elevation above the disc center h = 10. The 
observer is situated in the x-z plane at inclination i. For large inclinations and certain time delays, the isochrons consist of two disconnected regions: one is 
close to the flare, the other one is close to the disc edge (see e.g. panels b and c). At very large inclinations, a certain area on the disc surface becomes invisible 
(rightwards from the dashed curves given by eq. ka). Lower panels: Function t{r) = p — {z — h) cos i — r sin i. 



A2 Isochrons 

An a given moment of time, the observer receives the signal re- 
flected from the curve which is the interception of the paraboloid 
of equal delays with the disc surface. For an infinite flat disc, this 
curve is an ellipse with a center shifted along the x-axis: 

(a; — qsin^ i cosi)^ 



+ 



y 



(A3) 



(s/sini)2 (sctani)^ ' 

where q and s are given by equation (|3^). A few examples of 
isochrons for a flared disc are presented in Fig. |l^ At small in- 
clinations and small delays they are close to those for the flat disc 
(see Fig. |l^,c). The curve of equal delays on the disc surface can 
consist of two disconnected regions depending on the inclination, 
disc parameters and time. Some part of the disc can also be blocked 
from the observer by the disc outer edge. The area on the right side 
from the dashed curves on Fig. [l^ is invisible. 

A3 Integration over radius 

Computation of the response function for flared discs is reduced to 
one integral over the radius (see eqs. |l2| ^o|). For zero inclination 
i — 0, the integration is trivial using the (5-function from equa- 
tion (|2l|). For the non-zero inclination, the integration limits Vm, 
tm are given either by the solution of the equations 



or by the minimum/maximum disc radii rin, rout- Integrating over 
radius we use a substitution r — [tm + + (ru — fm) cos d]/2 
which removes possible divergency of the integrand at the integra- 
tion limits. The solution of equation (A4) is trivial in the case of 
z — (see eqs. ^0). In order to solve equation (A4) for 2 7^ 0, 
we first tabulate the function t{r) = p — {z — h) cosi — rsini 
at a dense grid of r, with r varying between —rout and rout- If 
tan i < Tout I H, equation ( A4 ) can have from 1 up to 3 solutions 
(for — r-out < r < Tout)- If t'(rovLt) > 0, then there are two solu- 
tions for iniin < t < t{r out) that are found by interpolation (tmin 
is the minimum value of t{r)). For rout) < t < t(rout), ru ~ 
Tout and the only solution is found by interpolation (Vm is then the 
absolute value of that). The response function for t > t(— rout) 
is zero. If t(rout) < (see Fig. [lo|b,c), the integral over radius 
has to be divided into two parts, since the isochrons consist of two 
separated regions for some t. 



For every t, we check that the visibility condition (^^2) is sat- 
isfied at all integration points 9. If the condition is not satisfied 
for some points, we find the radius (or radii) satisfying equation 
cos0(r) = cos 00 (r) by iterations (see eq s. [22 , A2). The integral 
is then recomputed only over radii where (K^T valid. All the rou- 
tines used in this paper are written on IDL. 



t = p — {z — h) cos i ± r sin i, 



(A4) 
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